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1 Introduction 



The differential equation 



1 1 1 

f{u) f{v) fiw) 
f{u) f{v) fiw) 



— 0, subject to u - 



(1.1) 



is one form of the first of a series of differential equations that can be written as 

n+l ^ 



1=1 



(1.2) 



subject to the constraint that / is an even function. These equations, which have appeared in the 
context of constructing mechanical systems with certain invariances, are the focus of this paper. 
(We will later review these connections between functional equations and mechanical systems.) 
When n — 2, (1.2) gives 

Odd 

if {Xl - X2) f {Xi - X3)) + — (/ {X2 - X3) / {X2 - Xi)) + 7^ (/ {Xs ~ Xi) f {X3 - X2)) = 0. 



dxi 



dx3 



(1.3) 



Using the evenness of / we may express this as 

d d 
- — (/ {xi ~ X2) f (x3 - xi)) + - — (/ (x2 - 2:3) / (a;i - X2)) 
ox I 0x2 



-Q^ if ix3 - Xl) f {X2 - X3)) = 0. 

(1.4) 



Equation (1.1) can be written in this form if we put u = xi — X2,v = X2 — X3 and w = X3 — xi, 
which automatically satisfies the constraint, although the assumption of evenness is not required. 

Braden and Byatt-Smith Q proved, as part of a more general theorem, that the complete 
solution set for the function / which satisfies (1.1) is 



/ (m) = a + bu or a + b e'^ 



(1.5a) 



bpicu + d, 32,53)- 



(1.5b) 



Here p is the Weierstrass p— function. Equation (1.5b) has six constants associated with it, 
namely {a, 6, c, d, §2, 53}, where 52 and 53 are the two constants which relate to the two periods of 
p (z), and d is one third of any period. However p (z) satisfies the equation p' = 4p^ — g2p — 33, 
with z^p(z) ^ 1 as z ^ 0. Hence p satisfies the scaling law p (cu, 52, .93) = c~^p {u,d^g2,c^g3) , 
so that without loss of generality we may take c = 1 in (1.5b). The solution set represented by 
(1.5a), which is a subset of (1.5b), does not require the constraint u + v + w = Otohe satisfied 
and is the general solution of the differential equation 



rr-r =o. 



(1.6) 



The solution set (1.5b) only satisfies (1.1) provided the constraint is satisfied and is the general 
solution of 



+ Af + Bf +Cf + D = 0, 
or, upon eliminating the arbitrary constants. A, B, C and Z?, 



(1.7) 



(1.8) 
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When we consider (1.5b) as the general solution of (1.7) or (pT^), d appears as an arbitrary 
constant. However if (1.5b) is also to satisfy ( pTl| ) then d is not arbitrary and by substituting the 
solution back into ( |l.l| ) can be shown to be either zero or any integer multiple of one third of any 
period of p{z). A simpler proof that all of the solutions of (1.1) are contained in the solution set 
of (1.5), subject to the above condition on d, can be obtained by eliminating the functions f (v) 
and / (w) by taking suitable combinations of derivatives of (1.1). At various points in the proof 
the equation factorises to give (1.6) or (1.7) as factors. (See Appendix A.) 

For the last 15 years the nature of the solutions to ( |1.2[ ) has remained ope n. One can show 
1^ that (1.5b), with c? = to ensure the evenness of /, satisfy equation (1^) and it has been 
conjectured that such were t he o nly solutions. Though the method of eliminating functions just 
noted above for the case of (1.1) should be applicable in the general case, the algebra involved 
to completely define the solution set is quite considerable. In the case of (1^) when n > 3 the 
amount of algebra appears to be so large that even a Maple calculation cannot handle the details 
necessary to provide a definition of the solution set. Here we shall develop new techniques to 
handle the equation. Our main result contains a surprise. To describe this let us introduce the 
following: 

Definition 1 Let A denote the class of meromorphic functions f{z) whose only singularity on 
the real axis is a double pole at the origin, and let Ap denote the class of periodic functions f{z) 
whose period is p (where p is real) and whose only singularities on the real axis are double poles 
at z — up, n G Z. 



Tiieorem 1 For functions f{z) G .4U Ap the general even solution of ( |Z.4 ) is: 

a) for all even n> 2 given by (1.5b) with d — while 

b) for odd n > 3 there are in addition to the solutions (1.5b) with d = the following: 



hi{z) = \/{p{z) -e2)(p(z) -es) = 



a^z) 



= b 



dn(w) 



9l{v) 4^26)3(0)614(0) sn2(M) 



h2iz) = V(p(z) -ei)(p(z) -63) 
hsiz) = Vipiz) - ei)(p(z) - 62) = 



f^iiz)(T3{z) 
aHz) 

ai{z)a2{z) _ 
a^iz) 



i{v)ei{v) 



9f(0) 



^ , cn(u) 

^liv) 4tj202(o)04(O) sn2(w)^ 



9f(0) 



H{v) 4c^202(O)(?3(O) 



cn(u) dn(u) 
sn2(-u) 



Here 



a{z+U}a) -zC(wc,) / 

(^a{z)^ -, °^ u=Vei-e3Z, 



63, 



with uj\ — Lo, UJ2 — —1-0 — Lo' and W3 — lo' , and we have given representations in terms of the 
Weierstrass elliptic functions, theta functions and the Jacobi elliptic functions vM. For appropriate 



ranges of z the solutions are real. These exhaust the even periodic solutions of (1.2) and their 
degenerations yield all the even solutions with only a double pole at a; = on the real axis. When 
n = 2,3 the theorem can be proved without the assumption that f{z) £ AiJ Ap provided it is 
assumed that / is meromorphic with a double pole at the origin. We also conjecture this latter 
assumption is all that is required for n > A but have been unable to prove this. The surprise 
is the appearance of these new solutions for odd n, which in turn yield new Poincare invariant 
mechanical systems. 



Our paper is arranged as follows. For completeness in §2 we will describe the origin of (1.2) and 



of the several connections between functional equations and mechanical systems. We then turn to 
methods for finding the solution set of (L2). In §3 two methods using a series solution introduced 
in Braden and Byatt-Smith |Q are discussed. These methods yield either a Laurent series for the 
solution set or a set of differential equations, whose common solution the solution set must satisfy. 
The advantage of these metho ds is that they define the solution set, or the differential equations 
that the set must satisfy, for (1.2) when n = 2 or ?i = 3. Even though these methods appear 
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intractable for larger values of n, the equalities derived for n — 2 and 3 turn out to be sufficient 
to define the solution set for all n. In order to establish this we require a new method, valid for 
all n. This is provided by the Fourier transform method. We develop this method for (1.1) in §4 
and §5 and for (1.2) for all n in §6 and §7. In so doing we prove theorem 1 and define completely 
these new set of solutions. 



2 Some Mechanical Systems 

Some years ago Ruijsenaars and Schneider |]l3| initiated the study of mechanical systems exhibiting 
an action of Poincare algebra 

{H,B} = P, {P,B} = H, {H,P} = 0. (2.1) 

Here H is the Hamiltonian of the system generating time-translations, P is a space-translation 
generator and B the generator of boosts. The models they discovered were found to posses other 
nice features: they were in fact integrable and a quantum version of them naturally existed. These 
models also appear in various field theoretic contexts (see Q). Ruijsenaars and Schneider began 
with the ansatz for a system of n -I- 1 particles interacting on the line, 

n+1 n+1 

H = coshpj ~ Xk), P — sinhpj JJ^ F{xj — Xk) 

j=l k^j 3 = 1 k^j 

and 

n+1 

B = -Y,x,. 

i=i 

With this ansatz and the canonical Poisson bracket {xi,pj} — Sij the first two Poisson brackets of 



(2.1) involving the boost operator B are automatically satisfied. Supposing further that F(x) 



F{—x), then the final Poisson bracket is equivalent to the functional equation 

n+1 

{i7,P} = 0^^a,n^'(2^j-2;fe)==0, (2.2) 

3 = 1 k^j 



With f{x) — F^[x) this is precisely (|l.2|), and so solutions of this equation (with even square 
root) yield Poincare invariant mechanical systems. At the time Ruijsenaars and Schneider were 
able to to show (1.5b) (with d — 0) gave solutions to these equations for all n. These solutions 
in fact yield n+1 independent, mutually Poisson commuting conserved quantities, and so are 
a completely integrable mechanical system. A scaling limit of the Ruijsenaars-Schneider model 
yields the Calogero-Moser system with Hamiltonian 

1 "^^ 1 

H = ^Y.P^ + qJ2p(i^-1i) (2.3) 

1=1 i=jL] 

which is another well-studied completely integrable system . 

We note that many connections exist between functional equations and integrable quantum and 
classical systems ^, ^, ^ 0, ^, ^, |ll|, |l^. Functional equation (1.1) (without any assumptions on 



the parity of the function /) arises for example when characterising quantum mechanical potentials 
whose ground state wavefunction (of a given form) is factorisable |l^. More recently it has 
been shown to characterise the Calogero-Moser system [|l[ . Several functional equations appearing 
in this setting and whose general solutions have still to be found are given in [||. 



3 



3 Series Solution Approaches 



We shall now describe two methods based on a series approximation for studying (1.2). 
Method 1: Obtaining a Series Solution. 

One method of attempting to prove the conjecture is to assume all the Xi,i — 1 ■ ■ ■ n + 1, are 
small and of the same order, so that we write write Xi — tQ. Then we assume that all the even 
functions, such as / (iC) can be expressed as a power series in t with C as an order one parameter, 
in the form 



2j-2 



(3.1) 



The constants Cj are given by = 1, if j < and Cj = 1/j!, if j > and are included for 
convenience. The series (5.1) allows for a double pole at the origin which can easily be shown 
to be the only allowable singularity. The coefhcients aj are determined by equating to zero the 
coefficients of the powers oft^^ in the subsequent expansion of (O). These coefficients are of course 



functions of Q as well as aj . However each coefficient factorises into a product of homogeneous 
polynomials in Q, independent of aj, and a factor dependent on the aj only. Equating this 
coefficient to zero successively determines aj for all j > 4 in terms of ao,ai,a2, and as which 
are arbitrary. This process can be completed to any desired order, J, if the expansion (^^) is 
truncated at a suitable finite value. Substitution of this finite polynomial into (1.7) shows for the 
cases where the method works, that for a suitable choice of {A, B. C, D}, (1.7) can be satisfied 
to any desired order. 

The proof of the above statement is, of course, incomplete: the form of the general term, aj, is 



not obtained and hence we cannot show that the full expansion ( |3T| ) satisfies {^). This method 
works for n = 2 and n = 4, but for values of n > 5 the amount of algebra involved becomes so 
large that even a Maple calculation cannot handle the details. When n — 3, this method does not 
work completely in that it leaves 04 arbitrary. Subsequent methods show that this is not the case 
and we will resolve this difficulty later. 

It is interesting to note that the same procedure works for solutions of (^]^), although the 
assumption of evenness is not required. Hence we look for a solution 



/(iC)=E« 



(to 



i-2 



(3.2) 



If we take oq 7^ the process automatically gives a2j+i = producing the same even function as 
obtained via (3.1) for the solutions of (1.2). However, if we take ao — then we find that ai = 
and (3.2) is then a Taylor series. We proceed as above to produce the coefficients Uj for all j >7 
in terms of 02 to uq and show that (1.7) can be satisfied to any order. However, 05 is not arbitrary 
and is given by 



03 05 



(3.3) 



This condition is automatically satisfied for all functions in the set (1.5a). However (3.3) is also 
equivalent to the condition p"~ {d) = 12p'^ (d) /p {d) , which is satisfied when d is one third or 
two thirds of any period of p (z) . This gives the required condition on the constant d appearing 
in (1.5b), when / belongs to this solution set. 

Method 2: Obtaining a Series of Differential Equations. 

An alternative method is to assume that one variable, for example xi is not small and write 
xi = X and Xi — tQ, 2 < i < n + I together with ( |3.lD , although xi can be replaced by any 
linear combination of the Xi, so long as n variables are scaled by t. The expansion of f{x — tC.) 
then naturally produces coefficients of P , which are function of / and its higher derivatives. The 



coefficients of P , in the expansion of (1.2) when equated to zero, now yield differential equations 



which must be satisfied by / (x) but contain the "arbitrary" constants, aj. 
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For the case n — 2 the first equation is 



aof" + 12ai/' - 12//' = 0. 



(3.4) 



This proves to be sufficient to determine a differential equation for / in the sense that the elim- 
ination of the constants uq and oi by differentiation gives (1.8). The higher order coefficients 



of P produce equations similar to (3.4) but contain the constants aj with j > 4 which are not 
arbitrary. In this case eliminating all the constants by differentiation will yield an equ ation which 
is still necessary but not sufficient to determine /. If we seek a series solution to (iA) by loolcing 
for a series solution of the form f{x) = X^j^o '^2^-22^^"'"^ then we obtain 60 — ao and bi — oi, 
with 62 and 63 arbitrary. The coefficients bj^j > 4 are then determined by these four constants, 
the recurrence relation being the same as in method 1. The equations which contain 02 and 03 
also require &2 = 02 and 63 = 03 and so also reproduce the series ( ^.l| ). 

This method is more co mplet e than the previous meth od i n that it does yield a necessary 
differential equation, namely (L8), that all even solutions of (|l.2|) for n = 2 must satisfy. However 
we cannot prove that all the differential equations produced by the expansion are satisfied. However 
this is not a problem since it is easy to verify, by substitution, that all even solutions of (1.5b) 
satisfy (1^) for this case. 



Again the method can be adapted to obtain the solution of (1.1) which are not even, using 



(3.2) with aQ — ai — 0. The procedure is the same except tliat because / (C) is not even we obtain 
differential equations containing / (C) and / (— C) • When one of these functions is eliminated we 
can then show that either (|1.6|) or (1.8) is satisfied. However, the differentiations required to 



eliminate the constants aj means that the necessary condition ( |3.3D is not recovered. 

For the case n — 3 the situation is more complicated since the series of equations similar to 



(3.4) only produce necessary equations for /. In order to show that / satisfies (1.7) we need to take 
two equations similar to (3.4) and find the set of solutions common to both differential equations. 
The algebraic details are quite complicated and require a Maple calculation. The details are not 
repeated here but a summary of the results is given. A copy of the Maple programme which 
produces these results with further explanation can be obtained by contacting the first author. 



The two equations, whose common solution satisfies (1.2), for n = 3 are given by (Al.l) and 
(Al.2) in appendix A. To determine the equation or equations for these common solutions we 
effectively eliminate the arbitrary constants ao,ai,a2 and a^. These constants are the ones that 
appear in the expansion ( |3.l|) used to derive (Al.l) and (Al.2). 

Before we derive these equations for the common so luti on we investigate the form of these 
solutions by looking for a common solution of the form (3T) but with the constants Uj replaced 
by bj. The result again shows that bo — aQ,bi — ai,b2 — 0,2 and — a^. However, initially 64 is 
not determined but 65 and 6g are determined in terms of bj,j < 4. However we get two different 
values of 67 and equating these values give an equation which has two solutions, namely 



63 = 



60 (262 _ ^^^2-) 



7bl 



or 04 



64 = 



60^2 



(3.5a,b) 



If we choose (3.5b) then in turn we find both equations require a common value for bj for j > 7 
and the series generates, as in the previous method, a solution of (1.7) or (1.8) and so belongs to 
the solution set (1.5). 

However, if we choose (3.5a) to be satisfied then the two equations have different solutions for 
bs and equating these values determines 64. This value differs from that given in (3.5b) and hence 
generates a solution / (z) = /i (z), with three arbitrary constants, which is not a solution of (1.7) 
or (1.8), unless 



66| 
56o 



(3.6) 



It is easy to verify that when (3.6) is satisfied, any finite truncation of the series is equal to the 
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same truncation of the solution 



/2(^) = foipfy^^,12,8) +61. (3.7) 



This is a subset of the general solution of (1^), which may be expressed as 





/2 (z) = 35i/ sin^ or h (^) - -36i/ sinh^ ' ' 

depending on whether is greater than, or less than, zero respectively. Since the choice of 

. °" . 

(3.5a) requires a specific choice of 63 and hence 03 it may be thought that this solution fails to 



be a solution of the full equation (1.2). However, as we will prove, this supposition proves to be 
erroneous. 

We can also determine the equation for the common solution of (Al.l) and (Al.2) by effectively 
eliminating the arbitrary constants ao,ai,a2 and 03. The details are contained in the Maple 
programme referred to earlier and are not given here. This results in equation (Al.3) which 
factorises, so that / must satisfy one of the following equations 

f — or f — constant, (3.8) 

clearly a solution, but not one of interest, and 

f"'f - f'f" = or /"// = constant. (3.9) 



We find that this equation only satisfies (Al.l) and (Al.2) if the constant appearing in ( |3.9[ ) is 
zero. Thus the only even solution again is / = constant. The third factor we recognise as equation 



(1.8) with solution (1.5b), with c? = for an even solution. The last factor is seen to be a fifth 
order ordinary differential equation which is cubic in /^"^. It is easily verified that the solution 
f ~ fi (z) defined earlier by taking the choice of (3.5a) for 63 satisfies this equation, (A1.4). 

At first sight it would seem unlikely that we could obtain the most general solution of (Al.4), 
which is even and has a double pole at z = 0. However, in §7 we show that the Fourier transform 
of the non-periodic solution is simple and easily inverted to give 

/ (z) = acoshfez/sinh^ 6z , (3.10) 

for arbitrary constants a and b. This coincides with the series for /i (z) obtained earlier as the 
solution of (Al.4) by choosing a = 661 and b — ^/%iJhQ, although this requires the choice 
62 = — 7a6^/60 = — 215^/56o- With this knowledge, one may guess that the periodic solution are 
of the form a cn (5z, k) /sv? {bz, k) a dn {bz, k) /sn^ (6z, k) or a cn {bz, k) dn [bz, k) /sv ? (bz, k), or 
the alternate forms given in theorem 1. These are all periodic equivalents of ( |3.10 ), where cn. 



dn and sn are the Jacobian elliptic functions and k the modulus. It is then possible to verify, by 
substituting into (Al.4) that the most general solution of this equation of the form required may 
be expressed as 

/i (z) = y b^W^ (z) + 26o6iI^ (z) -bl + ^5o&2, (3.11) 



where 



/ 20 {bob2 - 3bl) 861 (5bo62 - 3bf ) ^ 
Wiz) = p ^z, , j . (3.12) 

The quadratic in ( p.ll| ) has a double root when 66^ = 56o&2 which corresponds to the condition 



(3.6) which gives the function /2 (z) , defined by (p.7|). Thus /2 (z) is the two parameter family of 



common solution to ([Lq) and (Al.4), which includes the solution boz ^ in the limit bi 0. 



6 



If the equation for the p function is written as 



P 



4p^ - 52p - 33 = 4 (p - ei) (p - 62) 



63) , 



(3.13) 



in the usual notation, then the quadratic in (3.11) divides the cubic on the right hand side of 
(3.13). When ei and 63 are the common roots (3.11) can be simphfied to give the solution 
fi (z) — a en (bz, k) /sv? {bz, k) for suitable choices of a, 6 and k in terms of 601 ^1 s-nd 62- When 
62 and 63 are the common roots we recover the solution /i (z) = a dn {bz, k) /sv? (bz, k) while the 
third choice (ei, 62) gives /i = a cn {bz, k) dn [bz, k) /sv? {bz, k). 



4 The Fourier Transform Method 



We now look at the problem of determining the solution of (1.1) and ( |1.2| ) using Fourier transforms. 
There are a variety of difficulties which are not immediately apparent. These will be treated as 
they arise. They include requiring generalised Fourier transforms for functions of x which are 
unbounded either at infinity or at points on the real axis, and the consideration of distributional 



solutions both of ( |l.3| ), (1.4) and of the transformed equation. 

We first assume that the function / in ( |l.4| ) is not necessarily even but is bounded on the real 
axis. Then we define the double Fourier Transform F{k, I) = T[F]{k, I) of i^(a;i, a;2) by 



F{k,l)= 1 1 F{xi,X2)e~'''''^-'^^^dxidx2. (4.1) 



With 



d d 

F {X1,X2) = (/ {Xi - X2) f {X3 - Xi)) + — (/ {X2 - X3) / {xi - X2)) 



+ 1^ if {X3 - Xi) f {X2 ~ X3)) , (4.2) 

0x3 

we have 

F{k, I) = ifce-*('=+')^='7(-0 fi-l ik-^^'^+'^^V (-fc) f{k + I) 

-*(fc + Oe-^('=+')^^/(-fc)/(0, (4.3) 

where / (fc) = f {xi) e~^^''^dx\ is the Fourier transform of / (xi) . When / and /' are bounded 
on the real axis it is clear that all terms in ^"(2:1, 2:2) are also bounded for all {xi,X2) and hence 
F {xi,X2) = 0, for all such bounded solutions /. Thus the appropriate equation to determine / (fc) 
is F{k, I) — 0, which gives 

k?{~l) fi-k -l) + lf{k) f{k + l)-{k + l) f{-k) f{l) = 0. (4.4) 

However when / has singularities on the real axis F (a;i,a;2) is not defined at all points in the 
{xi,X2) plane and we find that F {xi,X2) — 0, everywhere except on a set of measure zero. In 
particular for the even solutions of (1.5b) this set of points is the lattice xi = xj, mod (A) and 
X2 = xz mod (A), where A is the real period of the function /. For such functions / we show, in 
§5, that F{xi,X2) ^ 0. We also show in §5 that this means that F = is not the appropriate 
equation and derive the correct equation. 



Equation (4.4) is a rather complicated functional equation when / (— fc) 7^ / (fc). A Taylor series 
method produces a three parameter family of solutions. Two parameters are as a consequence of 
the fact that if g(fc) is a solution so is / (fc) = ag{bk) for all constants a and b. This is as a result 
of the scaling symmetries of the original equation ( [4.4[ ) . So essentially there is a one parameter 
family of solutions. However, it is not easy to recognise the solution from its series. The method 
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which appears to give the most simple solution is the following. We decompose / into an even 
and odd function of the form 



/ (fc) = fi (k) + tkh (k) , 



(4.5) 



where /i and /2 are even functions of k and t can be ±1. Substituting this expression into ( |1.4| ) 
yields an equation of the form 



A {k, I) + tB {k, I) + i^C [k, 1)^0 
and since t can be either ±1 this gives two equations 

A + C = and B = 0. 



(4.6) 



(4.7) 



Now we assume that I is small and expand (4.7a,b) as a power series in /. Equating the 
coefficients of the powers of I to zero gives a series of differential equations involving /i {k) and 
/2 (k). ^'^From the first two equations we obtain 



and 



where 



2(/2(fc)-a2)/i(fc) = fcai/aA;) 
kaifi (fc) + ff (fc) + fcVI (fc) - fi (k) ai + 2ea2f2 (k) = 0, 



ai = /i (0) and a2 = /a (0) 



Eliminating fi{k) we obtain 



(2(/2-a2)/^'-/^') +4/2(/; 



802/2 



2ai) = 



and one integration yields 



3a?/^' = 4 (a2 - /2) {cal + fl + 802/2') , 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



where c is an arbitrary constant. This constant is the third parameter referred to above. Exami- 
nation of the cubic in ( 4.12| ) shows that for all (real) constants c other than c = or —4 we have 
an oscillatory solution for /2. These solutions must be rejected on the grounds that if the original 
function / (x) is bounded, / (k) must tend to zero as fc — > ± 00. When c = —4, /2 = 02 and 
is automatically satisfied and (|4.9|) then yields 



(4.13) 



/i(fc) = •\/3a2^cot (^■\/3a2k/ai 



Again this does not represent the Fourier transform of a bounded function / {x) . However, it 
does illustrate the general feature of all periodic solutions of (4.12). When c 7^ the solution of 
(4.12) with /2 (0) = a2 is an elliptic function which is even. Clearly, when c 7^ —4, /2(fc) — 02 
has zeros when A: = or an integer multiple of the period of f2- Since (4.12) implies that /2 = 
and /j when /2 =02, at each zero other than fc = 0, fi{k) has a simple pole. Hence /i has 
a periodic array of poles, apart from k = where the singularity is removable. This feature is 
illustrated by the function ap peari ng in ( 4.13| ). 

When c = we can solve ( 4.12 ) with /2 (0) = 02 to get 



/2(fc) 



802 



2cosh2aA: + 1' 



(4.14) 
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where a = 02/01, so that 



3a2k cosh ak 



sinhafc (2coshafc + 1) 
The two possibiUties for / (fc) are then 

7(fc) - /i {k) ± kh (fc) , 



giving 



/i(fc) 



602^6 



4afc 



- 1 



and /2 (fc) 



- 1 



(4.15) 



(4.16) 



(4.17) 



If we normalise by choosing 6a2 — nai and ai = —2, then /i is the Fourier transform of the 
function 1/ sinh^ {x — itt/S) and /2 is the Fourier transform of 1/ sinh^ {x + in/3) . Apart from the 
addition of a constant these are p functions, with periods (00, tt) , and a double pole at the origin. 

Thus, apart from the double pole at the origin and the shift by o ne t hird and two thirds of 
the imaginary period of the p function there is a unique solution of (|l.4| ) which is bounded on 
(—00, 00) and tends to zero at 00. An interesting limit is 02 — > which gives /2 = and fi — ai. 
This is also seen to be the only solution of (IS) and (O) when 02 = 0, subject to the condition 
that /i (0) — fli /2 (0) = 02 = and /2 ^ as fc = 00. The inverse of the Fourier transform 
/ (fc) = ai is no longer a function, but the distribution aid (x) which can be viewed as the limit 



aiS (x) 



lim 



0-2 



asinh (— x/a + vri/S) 



where naia = 602. 

The only other regular solutions of ( |1.4| ), on the real axis, are ones that are either not bounded 
at 00 or are oscillatory. Such functions have distributional Fourier transforms. For example, the 
Fourier transform of e*"^ is 2Tr6 (fc — a) and as a distribution / = 2tt6 [k — a) satisfies (4.4), since 



An'^k5 {-I -a) Si-k-l-a)+ Att^IS {k - a) 5 [k + I - a) 

-4:TT^{k + l)5{-k-a)6{l-a) = 0. (4.18) 

This is because k 6 {—I ~ a) 5 {—k — I — a) is zero unless l + a — l + k + a — Oovk — 0, I — —a, 
when the coefficient of the product of S functions is zero. A more formal proof can be obtained 
by defining the inner products {u,v)k,i and {u,v)k by 



{u,v) 



k,l 



u V dkdl 



and 



(m, = j u V dk, 



(4.19) 



with a similar definition for {u,v)i. Then 

{kd {-I -a)S{-k~l-a), F (fc, l))k,i = {S (-1 - a) 5 {-k - I - a) , kF (fc, l))k,i 

= {5{-l^a), -{l + a) F{-l~a,l))i (4.20) 

= 0. 



Although this only shows that e'"^ is a solution of (1.4) by analogy we can also have e"^ although 
this does not have a Fourier transform, except formally by analytic continuation. 



The Fourier transform of x is 2Tri6'{k) and using (4.19) it is straightforward to prove that 
kS' (-1) S' {-k - + ''5' (fc) 5' {k + I) and (fc + I) 5' (-fc) 5' (I) are identical distributions. For ex- 
ample 



{{k + l)5' {-k)5' {I), F{k,l))k,i = {S' {-k)S' {l),{k + l)F{k,l))k,i 

= {6' {-k),-F{k,0)-kFi{k,0))k 
= -/fc(0,0)-/z(0,0). 



(4.21) 
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A similar calculation shows that 



{kS' (~0 S' i~k F (fc, l))kA + {16' (k) 5'{k + l),F {k, l))k,i = -fk (0, 0) - / (0, 0) . (4.22) 



In addition we can also see that if / (fc) is not continuous / (0) can be arbitrary since (4.4) is 



satisfied automatically when either fc, I or k + l equals zero. In particular, we can add an arbitrary 



(1.4) 



multiple of S{k) to any solution of (4.4). This corresponds to adding a constant to any solution of 



The above distributional solutions cover the solutions of (1.5a) that do not satisfy (1.5b). To 
cater for the non even periodic solution of (1.5b) that are also bounded on the real axis, we assume 
that the solutions / (x) are 27r periodic and write 

P— — 00 

The factor 27r is introduced so that the Fourier Transform takes the form 

oo 

f{k)= apS{k-p). (4.24) 



Introducing this expression into (4^) and recognising that f{k) is only non zero when k is an 
integer, we obtain 

ka_LH~I' + L)a_K-L6{-k- 1 + K + L) + I {k - K) uk+l^ {k + I - K - L) 

~{k + l)a-K5{-k + K)aL5{l - L) = 0, (4.25) 

at all integers values of k and I. This is a distribution supported aX k = K and / = L and is 
identically zero if 

Ka_La_K-L + LaxaK+L - {K + L) a^KO-L = 0, (4-26) 
for all integer values of (K^L). Since fi{k) and /2 (fc) defined by ( 4.16| ) satisfy the continuous 



version of (4.26), it is clear that the solution ax — fi {K) , K ^ with gq arbitrary, will satisfy 
(4.26). It is also easy to construct the corresponding function / (x) that satisfies ( |l.4| ). If /i (x) is 
the function whose Fourier transform is fi (k) then we define 

oc 

ff(^)= E hi^-^^P)- (4-27) 

p— — OO 

^ oo 

This function is 27r-periodic and has Fourier series — ^ ax^-^^^ where 

aK^ E /i - 2?"^) ^''''"dx - / fix) e-'^'^dx - h (K) . (4.28) 

"'0 p=-oo 



The above solution for uk can also be obtained directly from ( 4.26 ) by successively solving all 
equations with 1 < \K\ < N, 1 < |L| < iV for = 1,2,3,..., the equations where K — or 
L = being automatically satisfied. Apart from the fact that oq is undetermined and is thus 
arbitrary, this process, as in the continuous case, produces a three parameter family of solutions 
with a_2ifi-i and ai arbitrary. If we choose a_i — ^o? j [o^ — l) and a\ — fia^ / (a^ — l) then 
the choice a„2 = /3a^/ (a^^ — l) yields aK = (iKa^^ / {a^^ — l) , which is clearly equivalent to 
/i {K) when [3 — 602 and a — exp {0,2/0,1) ■ Since all continuous solutions / {k) give a corresponding 



solution Ok — f {K), via ( [4.28[) , we presume that other choices of a_2 give solutions for ok which 
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are oscillatory and do not tend to zero as K ^ oo. As in the continuous case this gives Fourier 
series which do not come from a continuous function of x. 

It is also easy to show that there are additional solutions of (4.25). If the set Si = {o-k} solves 



(4.26) then so does the set S = {j^aK} provided = 1, for all integer values of K. This gives 
two additional solutions \e^'"^^l^aK\ and {e'^'^^^^^ax} ■ If /i {x) is the 27r-periodic function whose 
Fourier series is given by the set Si, then these two solutions are from the functions /i (a; + |7r) 
and fi (x + ^tt) that is the original function shifted by one third and two thirds of its period 
respectively. 

In this section we have shown, that by requiring / to be bounded on any compact subset of 
the real axis, we can recover all the solutions of (1.5a,b) which have this boundedness property, 
by taking the the appropriate Fourier Transform of (1.4). We have also shown that there are 
no other solutions / which have this boundedness property. However by taking the appropriate 
limit we have also shown that the 6 function is also a solution and have recovered this from the 
transformed equation. We have not obtained the p function solutions, which have double poles 
on the real axis. This is because the transforms of these functions do not satisfy the transformed 
equation (4.4). The resolution of this problem is the topic of the next section. 



5 Even Solutions 



kfii) + if{k)f{k + i)^{k + i)f{k)f{i). 



We now look at the solutions of (1.3) which correspond to the even solutions of (1.4). Equation 
(4.4) for the Fourier transform now reads 

(5.1) 

Clearly if / (fc) is a function not identically zero then ( |5.l| ) implies that k/ f (k) is linear so that 
/ (fc) is constant. This however only gives the distributional solution where / (x) is a constant 
multiple of S (x), although as a distribution / (fc) — a6 (fc) is also a solution so that / (fc) ~ aS {k)+b 
is a solution. These are the only even solutions obtained via this method. 

However, we know that the p functions, which have double poles at the origin, and in particular 
the function l/x^ satisfy (1.3). We do not recover these solutions from (5.1) because in these cases 
the left hand side of equation (1.3) is not identically zero but acts as a distribution. 

For example, if we substitute f = 1/x^ into the left hand side of ( |1.3| ) we get 

( \ 2 2 

g(Xl,X2,X3) = - 



(xi 


- X2)^ {Xl 

2 


- 2:3) 


(xi 


- X2 


f {xi ~ x^f 
2 




- 2:3)^ {X2 

2 


- Xl) 


{X2 


-xz 


f {X2 - Xlf 

2 


{X3 


- 2:1)^ {X3 


- Xl) 


{Xi 


-Xl 


f {X3 - X2f 



(5.2) 



This expression is identically zero outside the neighbourhood of xi~X2 — 0, xi — x^ = 0, X2—X3 — 
0. For fixed a;3 this is the region outside the lines xi — X2 = 0, xi = 3:3 and X2 — X3 . In Appendix 
B we show that (5.2) acts as a distribution supported only at the point xi = X2 = X3 in the 



{xi,X2) plane and can be represented by the distribution 

g {xi,X2,X3) = TT^ {Xi - X3) S" {X2 - X3) + S" {Xi - X3) S' {X2 - X3)^. 



(5.3) 



Hence if we take the Fourier transform of ( |l.4D assuming that / is even and has a double pole at 
the origin but is otherwise bounded, we have instead of = 0, 



F = 



7r^((5' {xi - X3) S" {x2 - X3) + 5" {xi - X3) 5' {x2 - X3) )xe 



-ikxi —ilx2 



dxidx2 



= iTT'^kl{k + l)e-''^''+^'^''\ 



(5.4) 
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Thus (4.4) becomes 

(fc/ (0 + if{k)) f{k + i)^{k + i) f{k) J {I) + ^^ki (k + i). 

The Fourier transform of is — tt |fc| and the identity 

{k\l\+l\k\)\k + l\ = {k + I) {\k\\l\ + kl) , 



(5.5) 



(5.6) 

for aU k, I ensures that the even functions / (fc) = ±7r |fc| satisfy ( |3.5[ ). We note here that if /(fc) 
satisfies ( ^.5|) so does — / (k). However, in addition to solving (^^) we also require that the inverse 
/ {x) satisfies the condition that x^f (x) ^ +1 as a; ^ 0. Thus although / (k) = ±7r \k\ satisfies 
( ^.5D w e must take the solution / (fc) — — vr |fc|. However, / (fc) = ±7r |fc| are not the only solutions 
ofp^. Writing (|5^) as 

7(0 (fc (7(fc + - /(fc)) - Ifik)) + lf{k) f{k + /) = n^kl (fc + I) (5.7) 
and then dividing by I and taking the limit as / — s- gives 

aokf (fc) - ao/ (fc) + ? (fc) = 7r2fc2, 
where ao — f (0) and the solution of ( ^.8D with / (0) = ao is 

/ (fc) = Trfc coth {nk/ao) . 
This is the Fourier transform of the function 

2 /aox^ 



(5.8) 



(5.9) 



/(^) 



1 



ao |ao| / sinh 



(5.10) 



Again, apart from the addition of a constant and a scaling of x this produces the unique p function 
of imaginary period tt and real period infinity. 

For functions / (x) which are 27r-periodic in addition to having a double pole of the form 
at the origin, we need to modify ( [3.3| ) and (5.4). The distribution g (xi,X2,X3) must be replaced 
by one which is repeated 2tt periodically in xi and X2 that is 



gi {xi,X2,X3) = ^Cp 



so that 



where 



47r2 CpqS {k-p)5{l- q) e 



-i{p+q)x3 



C, 



pq - -^PQip + 



Equation (5.5) now becomes 



(fc/(0 + ^/ (fc)) / (fc + - (fc + ./ (fc) 7(0 + E ^'^^ {K + l)5{k-K)5{l-L). 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



Again following ( 4.24 ) we write / (h) in the form 

oo 

7(fc)- apHk~p) 



(5.15) 



p—^oo 
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and obtain the following recurrence relation for the set {uk} 

[KaL + LaK) uk+l ^ {K + L) gkul + tt'^KL{K + L). (5.16) 



As in §4 with the recurrence relation ( 4.25| ), oq is arbitrary. Writing apc — TrKhx for K = 
1, 2, ... we obtain 

{b2 + bK)bK+L = hKhL + l (5.17) 

and with bi = — and L = 1 we have, 

1-/3 

bii+i = (6if (/3 + 1) + 1 - /3) / {bK {I -13) +(3 + 1). (5.18) 
This is easily solved to get 

= or OK = TT \K\ (^ ^_^|^| j . (5.19) 

This reproduces the result that gk — f [K) at integer values of K where / [K) is the continuous 
version. Also the Fourier series it \K\ tttitv | e'^^ can be recognised as the Fourier series 

-oo vi-r' '/ 

for the p functions again up to the addition of a constant. Here q = (3^^^ is the usual notation for 
the nome. 



The conclusion is that the only even solution of (1.4) are those of (1.5b) with d — 



6 The General Case 



We now consider the even solutions of (1.2) for a general integer n. It will be convenient to write 



(|1^ as 

g(x,Xn+i) = ^-^Y[f{xp-Xq), (6.1) 



„+l g 



q¥=p 



where / is even and x is the vector (xi,a;2,- • -Xn)- We then define the n-dimensional Fourier 
transform g(k, Xn+ij) = T[g\{\i, Xn+i) of g (x, Xn+i) by 

g{k,xn+i)= <7(x,x„+i)e-*'^''dx, (6.2) 

as a generalisation of ([4.l| ). Again there are problems with the double pole of / at the origin. 
Now however we find that g acts as a distribution over any plane through the origin and the 
singularities become more difficult to deal with. To overcome this difficulty we assume that the 
arguments of g are complex and we replace Xj by Xj + iej where e„+i > e„... > ei > and assume 
that e„+i is small. In other words 

ri+l g 

p=l P q=ip 



We then assume that in the definition of 5, (6.2), we integrate along the Real axis in the complex 
Xj +iyj plane. If the function / (z) has double poles on the x = Re z axis and no other singularities 
in the neighbourhood of the x-axis then provided Cn+i is small enough, there will be no singularities 
of g within the domain of integration of the integral occurring in (|6.2[) . 
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To calculate 5 we first write S^q = xi —Xn+i +« (ci — ^n+i) when q = 1, = xi— Xq + i {ei — e.q) 
when 2 < q < n + 1, and rjq — Xq + i (ei — e^), where Xq ^ xi — Xq. Noting that ^1 = ^„+i we 
consider 



^ n+l 

J J llf{^,)e-'''--dx"'-'dx. 



-1 9=2 



/(a)e 



—ikixi 



n / f{^,)e-''^''^dxq\dx. 



-ikiXi 



f n 

n 



(6.4) 



Now since > ei the singularities of / (zi), at zi = i(e„+i — ei) and / (zq), at = i(eg — ei), 
g > 1, lie in the upper half plane. 

For functions / (z) which are analytic in the neighbourhood of the Re z axis but have a pole 
at z = we define Fourier transforms fu and /l by 



fuik)= (pfix)e~'''^dx 



(6.5) 



and 



fL{k)= ^f{x)e-'^^dx, 



(6.6) 



which are suitably indented to go above and below the singularity at the origin. With these 
definitions we have, in the limit as e„+i 0, 



/(xi - x„+i + z (ei - e„+i)) e V'-i / da;i _[ _[ /l (~fcg 



q=2 



71 \ n 



\<;=1 / q=2 



(6.7) 



Thus the first term in the sum in ( 16.31) contributes a term ikil to g. A similar calculation gives 



^9=1 



V9=i+1 



V9=l 



,i=i / \9=i / 



so that 5 = gives 



E 



\9=1 / \q=J + l 



(6.8) 



/me^O- E^n n^^(-fc«)h'^- 



V9=l 



= 1 / \9=1 



(6.9) 
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We now define the Fourier transform of / (2) to be i (^fjj (k) + /l (k)^ . This corresponds to the 

usual generalised definition of the Fourier transform of functions with non-integrable singularities. 
This also coincides with the use in §5. Also we have the result that the difference Jl {k) — fjj (k) 
is 2wi multiplied by the residue of the functions / (z) e"*'^^ at the origin, assuming that the only 
singularity of / (z) on the real axis is at the origin. We also make the assumption that / (z) is an 
even function of z with a double pole of the form at the origin. Hence with 



and 



we have 



and 



/ (fc) = I {fu (k) + II (k) 



fh (k) - fu (k) = 2TTi X Residue (/ (z) e''^'-) = iirk 



h {k) = .f{k) + nk 



fu {k) ^f{k)~ TTk. 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



Hence the equation 5 = 0, (3.9) gives 



J = l 



9=1 



q=j + l 



\q=l 



9=1 



E^^- I [11 [fikq) + ^h)j = 0. (6.14) 



This equation is to be regarded as an equation which determines / (fc) and must be satisfied for 

all {kg} in M". Again we note that, by inspection, / (fc) — ±7r |fc| satisfies ( |6.14| ) and we may 

also show that f{k) = tTT\k\ satisfies (6.4) for the all values of t which satisfy (1 

jir 



tT:\k\ satisfies 
{t — 1)" {1 + t) . Apart from t = ±1, these are, t 



i cot ■ 



J 



1, 



<)"(t-l) = 
However, to 



n - 1' 

satisfy the requirement that the inverse / (z) is such that z^/ (z) ^ 1 as z — > we require the 
solution above with t = —1. This problem recurs throughout §6 and §7 and we will assume that 
we only take the multiple of / (k) which satisfies the criterion that it has the correct double pole 
either at the origin or at the sequence of double pol es w hen / (z) is periodic. With m = 2 it is a 
simple matter to verify that iS'2 = is equivalent to (^.5|). 



The following lemmas prove useful in finding the complete solution to ( 3.14 ) for all n. We 
begin with a definition. 



Definition 2 For each n let 
Tn be the solution set of Sn 
Qn be the solution set of Sn 



0, 

and f (0) 0. 



Lemma 1 For n > 4 and even, Tn ^ F^- 



Lemma 2 For n > 3, Gn Q G2- 

Proof of Lemma 1: In Sn+2 we put kn+2 



Sn 



+2 



-kn+i and using the fact that / is even we find 



= (?{kn+l)-ll^kl+^)Sn. 



(6.15) 
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The factor (fc) — Tr^fc^ produces the solutions / (fc) = ±7r \k\ if / is even. This function we have 
already shown (see (5.5), (5.6)) belongs to the solution set ^2- Since the solution set Tn+2 must be 

contained in the solution set of Sn+2 |fc„+2=i-fe„+i, ( Fl^ ) shows that J^„+2 C (J{k) = ± |fc|)uJ^„. 
The result now follows by induction. 

Proof of Lemma 2: In we put kn+i = and obtain 



n+l 



/(0)5„. 



(6.16) 



Again if /(O) ^ the result follows by induction. 



Definition 3 Let B denote the class of functions whose generalised Fourier transform f (k) arise 
from functions f{z) bounded on the real axis apart from a double pole of the form l/z^ at the 
origin. 



This is the class of functions for which we have derived ( |6.14 ) . We then have 
Theorem 2 For solutions with f (z B, then Tn ~ T2 for all even n. 

Proof: By Lemma 1 Tn Q ^2- In §5 we have proved that T2 is the one parameter family 
/ (fc) — Tik coth (7rfc/ao) , together with its limit as op 0, ±7r|fc|. It is easily verified by 
substitution that these satisfy Sn — 0. This gives T2 C JF„. The result now follows. 

Theorem 3 For solutions with f Cz B, then Q„ = Q2 for all n. 

The proof of theorem || follows by the same method as the proof of theorem |2| However we 
cannot exclude the possibility that there exist solutions of Sn = 0, with n > 3 and odd, which 
have / (0) = 0, but are not contained in the solutions of 52 = 0. Before we look at the possibility 
of such solutions we consider the extensions of Lemma 1 and 2 and the corresponding theorems 
about the solution sets, to functions / (z) which are periodic. This requires / (z) to have a periodic 
array of double poles of the form l/{z ~ 2mi) at the points z = 2mT, n e Z. Again we follow the 



method outlined in §4 and write / (fc) as in (4.24) and reproduce from 5^ = a similar equation 
with f [kq) replaced by au , which is equivalent to (4.26). The conclusions can be summed up 
in the following theorem, which incorporates the result of Braden and Byatt-Smith |Q for even 
functions as a special case. 

Definition 4 LetB2Tr denote the class of functions whose generalised Fourier transform f (k) arise 
from 27r periodic functions f{z) bounded on the real axis apart from double poles at l/{z — 27rn)^, 

n e z. 



Theorem 4 For f e BU 62^ 

1. The only even solutions of (l.i) with n even are those of (1.5b) with d = 0; 

2. The only even solutions of (l.i) with n odd and for which / (0) 7^ are those of (1.5b) with 
d = 0. 



7 New Solutions 



We now look at the solution of ( |l.2| ) which have / (0) = 0, to see if there are solutions which do 
not belong to the set J-2. Lemma 1 and Theorem ^ can easily be adapted to prove that when 
n is odd J^n Q ^3- We first find this solution set and then prove that when n is odd ^„ — Tj,. 
Hence we look for solutions of (6.14), with / (0) = 0, when n = 3. We wish to consider only even 
functions, but wish to include functions like |fc| which is not differentiable at fc = 0. Hence we 



consider (6.14) defined on the subspace fcg > 0, g = 1,2 and 3, with all derivatives at the origin 
defined by one sided derivatives. Thus in the interval fc > 0, |fc| is defined as k with derivative 1 
at the origin. 
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In (6.14) we write ki — k and k2 = k-s = I and let / 0. Then (6.14) gives 
flo (aokf (k) - aofik) + ? (fc) - ir^k^) = 0, 



(7.1) 



where uq = f (0) . When ao ^ this gives the same equation as (^^) but is automaticaUy satisfied 
when oo = 0. So in addition to the solution given in (5.9), which belongs to J^2, we can also allow 
flo = 0. 



When oo = the next term in the expansion of (6.14) gives 

ai (? (fc)-^'fc2) =0, 



(7.2) 



where ai = /' (0). If oi 7^ then (7.2) gives / (fc) ~ ±7r |fc| as the only even solution. This is also 
the solution of ( |5.9| ) when — and also belongs to ^2- Now if we assume that ai — 0, we can 
write the third term in the expansion of (S.14) as 



2Tr^kf' (fc) + 02/^ (fc) - 27r2/ (fc) = 7r^a2k^, 
where 02 = /" (0) . The only even solution of this equation is 

f (k) = irk tanh [ 

\ 27r 



(7.3) 



(7.4) 



which automatically has /" (0) — 02 . This of course is a necessary requirement and we need to 
check that this is a s olution of (|6.14 ). 
We rewrite ( |6.14 ) as 



/ [kg) - Trfcg 

{ f (kj) + TTfcj ^^^^ \fikq) + nkq^ 



n 



f\EkA+7Tj:k, 



-1 = 0. 



Substituting / (fc) = nk tanh (ak) into 5„ gives 



(7.5) 



Sn-ij: (e-'^^ + 1) n 



1=0 



The first term can be written as 



exp 2a ^ kj 

/ „ \ 
exp 2a X] 



S^n'' = ^ ftj+i - flj = a„+i - ai, 



where = (-e-^^fc-^) so that S^i^^ = (-1)"+' exp I -2a ^ /cj + 1. Hence 



n+l 



1. 



9=J 



- 1 



exp 2a ^ fcj + 1 



(7.6) 



(7.7) 



(7.8) 



This gives immediately Sn = whenever n is odd. Hence / (fc) = tt/c tanh (afc) , with a arbitrary, 
is a solution for all equations Sn = 0, when n is odd. It is also evident from ( |7.8| ) that this solution 
does not satisfy S'„ = for n even. 
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We also note that substituting f (k) — 7Tkcoth{ak) into (6.5) changes (7.6) to 



Sn 




2akj 



exp {2akj) 

-2afe„ I 

1 / „ \ 

exp J2 "^-akj — 1 



(7.9) 



The change in sig ns no w means that 5„ = for all n showing that fi{k) — nkcoth (ak), a 
arbitrary satisfies ( 6.14 ) for all n as indicated earlier. 

If we write a = ^n/a so that / (k) = Trfctanh (^^irk/aj then this is the Fourier transform of 
the function f (z) = — ajajcosh az/sinh^ az, so a must be negative to satisfy the pole condition 
that z^f (z) —> 1 as z ^ 0. However, since the coefficient of the double pole is in fact arbitrary we 
have / (z) ~ /3coshaz/sinh^ az satisfies (1.2) for all odd values of n. This is the solution referred 
to in §2. 

The even period solutions, which satisfy the modification of ( |6.14 ) when / (z) has an array of 
double poles of the form (z — 2p7r) ^ at the points z = pir, p = 0, ±1, ±2, can be written as 



/ (fc) = QpS {k — p) , with Q-p — Up, 



(7.10) 



P— — 00 



as in ( 4.24 ). Again, ( ^.14| ) is now to be satisfied at all integer values of {kg} with / (kg) replaced 
by flfc . Hence 



9=1 



TTNg) II {aN,~7rNg) 



71 



71 



Y[{aN,+7rNg)\ ^0. (7.11) 



To solve 5*3 = we proceed as in §5. If ag is not equal to zero we can put — and recover the 
solutions ( ^.19 ). However, if op = then 5 = for all Ni and N2 if A^3 = 0. Writing down all 
the equations for Nj > 1 it is easy to establish that if ai = tt {(3 + 1) / (1 — /3) the odd terms are 
given by 



a2N- 



_,=7r(2jV+l) \_^2nJ ' ^>0' 



(7.12) 



as in (5.19) but the even terms depend on the choice of 02 which must take one of the values 

(7.13) 



(I + /32) (1-/3') 
a2 = 27r-!^ -J- or 27r-^ ^ 



1-/32 



1 + /32 



The first choice gives 



and the second 



(l + /32^) 



1-/3: 



(1 - p^^) 
a2N^2TrN\ iV > 0. 



1 + /3 



2N 



(7.14) 



(7.15) 
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The result (7.13a) and (7.14) is equivalent to (5.19) while writing /3 
is equivalent to 



aN=TT\N\ 



l-/3l^l 
1 + /3I 



-/3 (7.13b) and ( |7l5| ) 



for all N. 



(7.16) 



This reproduces the result of (7^) at integer values of k. T he proof that {ajv} satisfies 5„ for all 
n is identical to the proof in the continuous case, (see (|7.5| ) - ([z!^)). 

The corresponding inverse / (z), obtained from /(fc) can be constructed either by taking the 
Fourier inverse of / (fc) or by the infinite sum defined by (3.27) using the function /3 coshaz/ sinh^ az. 
As noted in §3, this function must be one of cn/sn^, dn/sn^ or cndn/sn^. There arc two reasons 
why there are three functions representing the solution set. The first is that if all the parameters 
defining the elliptic function are real then the transformation z iz permutes these three func- 
tions according to Jacobi's Imaginary transformation sn(iz, k) i sn(z, 1 — /c^) /cn(z, 1 — /c^) , 
cn{iz, k) — > l/cn[^z, 1 — fc^) and dn(iz, k) ^dn(z, 1 — fc^) /cn(z, 1 — fc^) . Secondly, Jacobi's Real 
transformation defines the elliptic function for the parameter fc < or fc > 1 in terms of elliptic 
functions with a scaled independent variable and parameter fc in the range < fc < 1. Again the 
effect is to permute the three functions. 

The conclusion is that the even solutions of (]1.4|) with n odd fall into two categories. The set 
defined by (1.5b) with d = 0, which satisfy (|1.8| ) subject to z^f (z) a constant as z ^ and 
also the set which satisfy (A1.4) subject to the same condition at the origin. These solutions have 
been expressed in terms of the p function and also in terms of the Jacobian elliptic functions and 
theorem 1 is proved. 
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Appendix A 

The expansion procedure, for 71 = 3, described in §3 yields the following two simultaneous equations 
(Al.l) and (A1.2) for the function f (x), which satisfies (|1.2|) 



120 02 /' / + 5 ao /"' /" + 60 /" /' ai - ao f^"^ / = (Al.l) 

and 

504 as f f + 1080 /' as /" / + 36 /' ai / + 15 /' aq /'"^ + 180 /' /"^ oi + 180 /"' /'^ ai - 
3 ao /''+360 /'^ a2 + 15 /' /^^^ /"-60 /" ai /"' /~12 ao /(™) /"' /+240 /"' P+ao /(^) /" / = 
0. (A1.2) 

These two equations determine / in the sense that the complete set of solutions for (^^) is 
contained in the complete set of solutions which satisfy both (Al.l) and (Al.2). The common set 
of solution to (Al.l) and (Al.2) satisfy the following equation 

jyf,, j„2_^,„2 ^, J,,, j(iv)^j{v) j,,^^ /"')(2000 p f'^ /" V125O f f'^ f'"^- 

18500 p f'^ /"^+50 p f'"^ /('')-250 /'"^ /" /s+63000 / /'^-3000 p f'"^ /"V18 p p"^^ 
70875 p /'"^ + 22500 p f'^ p'^ + 1692 p f'^ /(™)^ + 150 p /"' /' p"'^^ p"^ - 70875 /"^ + 
930 p f'^f^^^ /(™)^/"— 5220 /'"^ P^ f'^"^ f" /(™'+2475 j/4 j(t,) y"_33oo j'"^ j/3 j„2 

32250 /3 /'^ /"^-10125 /'^ /"^ /(")-27000 /'^ /"^ /(™)-2875 /'^ /" /"'''-9875 p /'^ f"^- 
174375 p p^ p'^ /'"+212625 / /'^° /"^ /'"+1875 f'^ /"^ /'"+182250 /'^ /"^ /"' ^31000 /'^ /"'- 
1500 /5 /"^ /'^ /"'-93375 p f'^ f" /"'^-8000 p f'^ /"^ /"' ^4750 /'^ /"^ /"'^+5625 /' /"^ f'^- 
2875 /5 /'"^ /"^ /"'^+375 /' /"^ /"''^-212625 p f" /"'^-60 p f"^ f /(")%90 /("^^ /"^ /"'- 
300 /' /^™' /'"''+750 /'"^ /"^ /f^^^+lSO /"^ /"'^— 900 p f" f'^ /(™)^— 144 j'^ /('")^— 

750 /"^ /'"^ /'^ /(")+5190 /'' /'"^ /'^ /(™)^+90 /"^ /"' /(")^-1215 /'' /'^ /"-441 /"^ /(™)^ /*^'"^ + 
30 p f'^ /(")^ /"^ — 14250 p f'"^ f'^ 90 /'^ /" /(^)^+1875 /'^ f"^ /(™)^ — 10125 /'"^ /'^ Z^"^ — 

4185 p /"' Z'*^ /("')^-1300 p f'^ p^ /(")+1500 p /'^ /(™) /"^+2355 /'^ /(™)^ /"^-27000 /'"^ /'^ /(™) + 
4185 /''' /(™)^ /" — 9075 /'^ /"'' + 6375 p f'^ f'-^^ f"^ — 1815 /"' /'^ /" /(™) /'"^ + 
3240 /"^ /"' /'^ /"^ /^™)^+750 /'' /"' f""^ f /("^ — 750 /'' /'^ /^™) /"^ /(")+5220 /"' f'^ /(™) + 

300 /''' /'"'^ /''^ /"'^ /^^^ 7500 f'" f''^ f"'^ /^^^^+25 p /'"^ /^^'^ /""^ 12180 Z*^ f'^ p'^^)'^ /"~t~ 

2865 p P^ f^^^ /"^ /(™) + 18600 /"' /(™^ /"^ 8100 /"' /^^'^ /'^ /"^+120 f'"^ f'^ f'^^^ p^-") + 
150 /'"^ /' /(™)^ 1080 p f" /(^' /"^ /' /(™) + 120 p /"' /("^^ /'^ /"— 150 p p"^ f" /("^ — 
90 /9 /' /(^)^ /" /(™'+54000 p /"' /" /("') +20250 p /"' /'"^ /"+63900 /'"^ /'^ /"+ 
250 /'"'^ /^"' /' /"— 81900 /'™^ /"^ /"'— 1875 /'"'^ /'^ /" /'™)+975 /'"^ /' f"^ /(™)) = 
0. (A1.3) 

Equation (A. 13) has several factors. The only factor which provides solutions to (1.2), whenn = 
3 which are not of the form (1.5b) is 

(2000 p f'^ /" V125O p /'^ /"'^-18500 p /'^ /"^-250 /'"'^ /" 7^+63000 / p^ /'^-3000 p f"^ /"'*+ 
70875 p /"'^+22500 p /'^ /"'Vl692 p f'^ /(*")^-70875 /"^ /'"-3300 p /'"^ /'^ /"^ /(™) + 
32250 /3 /'■^ /(™) /"^-27000 p f'^ /"^ /(™)-2875 p f'^ /" /"''*-9875 /'^ /""^ /"'^-174375 /"^ /"'+ 
212625 / /'^° /"^ /'"+1875 /'^ /"^ /'"+182250 p /''^ /"^ /"' ^31000 p f'^ f"'^ /'"-1500 p f"^ f'^ /"'- 
93375 p /'^ /" /"'^-8000 p f'^ /"^ /"' ^4750 p f'^ /"^ /"'^+5625 p /' /"^ /"'^-2875 /''^ /"^ 
375 p p /"^ /"''*-212625 /'^ /" /"'^-300 /' /'"''+750 /'"^ /"^ /(*")-900 p f" /'^ /(™)^+ 
5190 p f'"^ p^ /f^")^— 1215 /'' /'^ /"— 14250 f'"^ f'^ /(™'+1875 /'^ f"^ /(™'^— 4185 /"' /'^ /(™)^+ 
1500 /'^ /(™) /"^+2355 p f'^ /(™)^ /"^-27000 p /'"^ /'^ /(''')+4185 J'' f'^ p'"^^ /"-9075 /""^ p^ + 
3240 /"^ /"' /'^ /"^ /("')^~7500 /"' /'^ /(™)-12180 /"' f'^ J^''')^ /"+18600 /"' /(™) /"^ /'**+ 
150 /'"^ /' /(™)^ /"+54000 /"' /'^ /" /(*^) +63900 f'"^ f'^ /^™^ /"— 81900 Z'' /'^ J^™^ /"^ /"'— 
1875 p p'^ p^ /" /("')+975 p /'"^ /' /"^ /("')) + (-750 p f'"^ /'^-1300 /"'^-10125 p /'^ /"V 
150 p p'^ p'^ -10125 p /'"' /'^-441 p f'^ /(^^■)'+6375 /'^ /"^+50 /"'Vl50 /"' /' /(")' + 
930 /8 /'^ /(^"^^ /"-5220 /5 f*^ /" /(*^)+2475 /'"^ J''' /"+750 p /"' /'-750 p /'^ /(™) /"^+ 
5220 p f" f'^ /(™)+300 /'"^ /'^ /" V25 /"' p^ /"^+2865 /"^ /'^ /"^ /(™)-8100 p f" f'^ /"V 
120 p /'"^ /'^ /(™)-150 p /'"^ /" /(*^') +20250 /"' /''^ /"+250 /'"^ /' /"-1815 p /"' /'^ /" /(™)- 
1080 /"' /"^ /' /(*"))/(■") + (—90 /' /" Z'^™^ + 30 p f'^ f"^ — 144 p f'^ /'"'^ — 90 p f'^ f" + 

90 P /"' - 60 p /' + 120 p /"' /" ^§0 p /"' /'")/(")' + 

IS f9 f(")^ = n ( ^^ A\ 



Appendix B 



In this appendix we show that g (xi, X2, x^) given by (5.2) acts as a distribution. For convenience 
we write Xi — x, X2 = y and x^ ~ 0, so that we consider 



where 



9 (x, y) = 91 {x, y) + 92 [x, y) + gz{x, y), 

— 9 9 

91 {x,y) = 

92 {x, y) 



{x — y) x^ {x — y) x^ 
-2 2 



and 



We wish to prove that 



53 {x, y) 



y3 {x - y) y^ {x - yf 
2 2 



2-3^2 x'^y^ 



g{x,y)^7r' {S" (x) S' {y) + S" (y) S' (x)) . 



(Bl.l) 
(B1.2) 
(B1.3) 

(B1.4) 
(B1.5) 



To do this we first regularise integrals which have singularities, using the Hademard finite part 
interpretation of these integrals. We will see that this is consistent with the usual definition of 
generalised Fourier transforms used in §4 and also with the definition of / (fc) in terms of fjj (fc) 
and /l (k) , see ( |6.10 ). 

For a function (f) (x) , which is bounded in the neighbourhood of the origin, we define 



dx = lim 

C-+0 



and 



^ i-oo 



<t>{x) 



dx — lim 



^dx - 



(B1.6) 



When (f) (x) = e we have, from the above definitions 



dx = — TT I A: I 



corresponding to the interpretation of the Fourier transform of l/x^ in §3 and §5. 
Now we let (j){x,y) be a test function and define 



(B1.7) 



(Bl. 



x^ {x — y)' 



rdy = lim 



£-►0 



4>{x,y) 2(t)2{x,x) 

^^2/ + 7— ;r 



\x-y\> 



x^ [x - y) 



e X'- 



(B1.9) 



where (\)2 {x,y) — t^, and hence 
oy 



PQO POO 



00 — 00 



— '^^ ' ^'^ 3 dxdy = lim | 
x'^ {x — y) 



a;|>e|,x — y|>e 

0(0^, ^2 

3—02/ + - 

yi e 



a;2 (a; — y) 



(Bl.lO) 



\y\><^ 



\x\>t 
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Similarly we define 



(f>ix,y) f (j){x,y) 2 4>{x,x) 



i™i / '"".^ dy--^^^]- (Bi.ii) 



Now since 



x^ {x — y) I J x^ {x — y)" ' e X' 

\ix-y\>e 



6{x,y)dy f (j){x,ri + x)dr] _ 1 mi 1o^ 

2= / 372 = 73^' 



a;3 {x - y) J x-^V 

\x,y\>€ |r)|>e 



5/ _ /■ {<f)i {x,r] + x) + (f)2{x,'n + x))dr] 



with 



dx J ' 
l'?l>« 

a 

where again {x, y) = —0 (a;, y), we can, with the use of (B1.17) define 



(B1.13) 



-oo J-oo (a; - y) J J (a; - ?/) 

|x|>e|x-y|,>e 



dxdy 



_2_ J J + 1(^,(0,0) + 0,(0,0)). (B1.14) 

|r,|>6 |x|>e 

Using similar definitions for the remaining four integrals arising in JJ (pd we may define 
J J (pgdxdy = lim ( J J (f>gidxdy + J J <f>g2dxdy + j j (j)gsdxdy 

M2 ^ Hi 02 fis 

+ 1 (</.(e, e) - 0(-e, -e) - 0) + </>(-£, 0) - 0(0, e) + 0(0, -e)) j (B1.15) 

This requires a good deal of simplification in integration by parts. The regions Q,\,Q.2 and Q.^ 
appearing in (B1.15) are given by 

fii= {|a;|>e , \x-y\>e , |y| < e} , 

^2= {\x-y\>e , \y\>e , \x\ < e] , (B1.16) 

= {kl > e , |y| > e , |a; - y| < e} . 

This simplifies to 

A lengthy, but straightforward calculation shows that the limit of the sum is 

(77^-4) {(l>^^y (0,0) + 0j,^^ (0,0)). 

Hence 

llcljg^T:^ {(j),,y (0, 0) + <i>yy^ (0, 0)) . (B1.18) 
This confirms that the function g {x,y) acts as the distribution given by (B1.5). 
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